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ABSTRACT 


/ 


The  construction  and  operation  of  a  chemical  generator  of 
O^(a^Ag)  are  described.  This  system  could  be  readil^  modified  to  drive 
a  purely  chemical  iodine  laser  operating  at  1.315  pm. 


Optical  techniques  were  developed  for  estimating  the  excited 
oxygen  concentration  that  arrives  at  a  downstream  cavity  at  total  0^ 
pressures  from  1  to  6  torr. 

Mechanisms  that  explain  the  deactivation  of  excited  0'2  in  the 
gas  stream  are  discussed.  A  value  for  the  singlet  delta  energy-pooling 
rate  constant  is  determined  from  the  measured  02(a*Ag)  concentrations. 


RESUME 


On  d&crlt  la  construction  et  le  fonctionnement  d'un  g£n£rateur 
chlmlque  d'oxygdne  02(a*Ag).  Ce  systSme  pourrait  Stre  facilement  modi- 
fi 6  afin  d'actionner  un  laser  i  lode  entllrement  chlmlque  qui  fonc- 
tionne  3  1.315  pm. 

On  a  d£velopp€  des  techniques  optlques  afin  d'estimer  la  con¬ 
centration  de  l'oxygdne  excitfi  qui  arrive  dans  une  cavi t§  situge  en 
aval  du  g£n£rateur.  Des  presslons  totales  d'oxygdne  de  1  &  6  torrs  ont 
£tudl(es. 

On  dlscute  des  m&canlsmes  qui  expllquent  la  dSsactlvation  de 
1' oxyg^ne  exciti  dans  le  flot  de  gaz.  A  partir  des  concentrations 
mesurSes  de  02(a*Ag),  on  determine  une  constante  de  vitesse  de  d§sacti- 
vatlon  par  un  processus  appeli  mise  en  comoun  de  l'§nergie  (energy- 
pooling)  . 
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1.0  INTRODUCTION 


The  three  lowest  bound  electronic  states  of  oxygen  are  denoted, 
in  order  of  increasing  energy,  02(X3£g),  O^a^Ag)  and  02(bl£g).  For 
brevity  these  will  be  referred  to  as  02(X),  02(a)  and  02(b)  respec¬ 
tively  throughout  this  report. 

The  first  excited  state,  often  called  singlet  delta  oxygen,  may 
be  generated  in  a  microwave  discharge  of  pure  oxygen,  or  by  bubbling 
chlorine  gas  through  an  aqueous  solution  of  basic  hydrogen  peroxide 
(for  a  review  see  Ref.  1).  In  1978  it  was  first  demonstrated  (Ref.  2) 
that  such  chemically  generated  02(a)  can  provide  the  energy  required  to 
sustain  iodine  atom  lasing  at  1315  nm.  The  potential  of  this  purely 
chemical  oxygen  iodine  laser  (COIL)  as  a  highly  efficient,  high-power 
source  of  directed  energy  has  led  to  several  reports  of  Improvements  in 
both  the  oxygen  generator  and  laser  output  levels  (Refs.  3-7),  and  to 
an  assessment  of  its  capability  to  drive  the  nuclear  fusion  process 
(Refs.  8,  9). 

The  fundamental  process  in  the  COIL  is  an  exchange  of  electronic 
energy  (E-E  transfer)  during  collisions  between  02(a)  and  ground  state 
iodine  I(2P3/2)»  producing  02(X)  and  excited  I( ZP 1/2) •  This  E-E  trans¬ 
fer  is  fast  (Ref.  10)  because  there  is  a  near  resonance  between  the 
participating  electronic  levels  of  the  two  species,  and  it  is  this 
speed  that  makes  the  population  inversion  in  iodine  possible.  Propor¬ 
tions  of  02(a)  in  excess  of  about  20%  of  the  total  oxygen  are  required 
to  operate  the  iodine  laser  above  threshold,  and  water  vapour  must  be 
minimized  as  it  is  an  efficient  quencher.  This  threshold  minimum  pre¬ 
cludes  the  use  of  the  ndcrcwave  discharge  as  an  02(a)  generator  in  the 
I  laser,  since  proportions  of  only  5-15%  are  attainable. 

High  energy  storage  densities  are  possible  in  02(a).  Firstly, 
because  the  chemical  generation  method  is  remarkably  efficient,  both  in 
terms  of  reactant  conversion  and  product  purity,  it  gives  close  to 
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100%  02(a)  at  the  generator.  Secondly,  the  gas  can  be  efficiently 
transported  from  generator  to  laser  cavity  in  a  vacuum  flow  system* 
02(a)  has  a  very  long  radiative  lifetime  (Ref.  11),  and  is  not  rapidly 
quenched  ty  itself  (energy-pooling) ,  water  vapour,  or  inert  carrier 
gases  (e.g.  argon). 

Downstream  02(a)  percentages  have  been  reported  (Refs.  3-8)  in 
the  25-60%  range  at  pressures  in  the  region  of  1  torr  total  oxygen. 

The  proportion  of  02(a)  arriving  at  a  measurement  cavity  depends  on 
several  physical  factors,  principally:  the  pumping  speed,  l.e.  the 
transport  time  from  the  generator;  the  total  02  pressure;  the  condition 
and  material  of  the  walls;  the  temperature  of  any  trap  which  may  be 
installed  to  remove  R^O  vapour;  and  the  presence  of  other  gases  in  the 
system. 


Work  was  initiated  at  DREV  to  develop  locally  the  technology  for 
the  chemical  generation  of  02(a)  in  order  to  provide  a)  the  capability 
of  building  an  iodine  laser  should  the  need  arise,  and  b)  an  energy 
transfer  medium  for  other  potential  gas  laser  systems.  In  particular 
we  studied  the  energy  transfer  to  gaseous  nitrogen  fluoride.  The  elec¬ 
tronically  excited  NF  so  produced  emits  at  529  nm.  This  application  of 
the  chemical  02(a)  generator  to  a  potential  blue-green  laser  gain 
medium  will  be  the  subject  of  future  reports. 

This  report  describes  the  operation  of  a  prototype  chemical 
generator  of  oxy gen  singlet  delta  (Chapter  2.0),  and  the  measurement  of 
02(a)  concentration  at  a  downstream  cavity  (Chapter  3.0)  by  optical 
methods.  In  Chapter  4.0  the  deactivation  processes  are  discussed,  and 
a  value  for  the  energy- poo ling  rate  constant  of  02(a)  is  calculated 
from  the  experimental  data  of  Chapter  3.0.  The  energy  pooling  rate  is 
Important  in  modelling  the  maximun  energy  storage  capability  of  singlet 
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delta  oxygen.  Its  value  Influences  chemical  oxygen  Iodine  laser  output 
and  efficiency  predictions,  and  hence  also  the  conclusions  concerning 
its  candidacy  as  a  fusion  driver. 

This  work  was  performed  at  DREV  between  January  1981  and 
October  1982  under  PCN  33H07,  Research  on  Chemically  Excited  Lasers. 

2.0  CHEMICAL  GENERATION 


This  chapter  describes  the  construction  and  operation  of  a  chem¬ 
ical  generator  of  02(a)  similar  to  those  described  in  Refs.  4  and  5. 

The  excited  oxygen  is  generated  by  bubbling  chlorine  gas  into  an 
alkaline  aqueous  solution  of  hydrogen  peroxide.  The  overall  reaction 
is: 


Ci2  +  20H“  +  H202  -  2CJT  +  2H20  +  02(a)  [l] 

The  mechanism  of  this  solution  reaction  is  discussed  in  the 
literature  (Refs.  6,  12,  13),  and  will  not  be  considered  here.  It  is 
sufficient  to  know  that  it  is  a  very  fast,  essentially  stoichiometric 
conversion,  and  that  virtually  no  chlorine  leaves  the  solution  under 
the  flow  conditions  quoted  here. 

2.1  Apparatus  and  Operation 

Figure  1  is  a  schematic  representation  of  the  flew  system  which 
we  have  used  to  generate  02(a). 

Vessels  C  to  F  and  H  were  of  Pyrex,  as  was  the  tubing  between  D 
and  G.  Where  possible,  tubing  was  of  large  diameter  (40  mm)  to  mini¬ 
mise  02(a)  losses  by  wall  collision.  The  measurement  cavity  G  was  of 
rectangular  cross  section  (1  x  12  cm),  and  about  30  cm  long.  It  was 
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F  Hater  vapour  trap 
G  Measurement  cavity 
H  Liquid  nitrogen  trap 
I  Vacuum  pump 
T  Thermometer 

FIGURE  I  -  Flaw  system  for  generation  of  oxygen  singlet  delta 


made  of  Perspex  coated  on  the  Inside  with  Teflon.  Plastic  was  used, 
and  the  cavity  was  fitted  with  aluminum  electrodes,  so  that  high  volt¬ 
age  discharges  could  be  applied  transverse  to  gas  flows  containing 
02(a),  nitrogen  fluorides,  and  other  gases.  Gas  flow  rates  and  pres¬ 
sures  were  measured  with  calibrated  Hastings  flowmeters  and  capacitance 
manometers  (MRS  Baratron,  Vacuum  General)  respectively.  The  vacuum 
pump  was  an  Edwards  660,  which  is  specified  as  being  capable  of  dis¬ 
placing  11  L/s. 

A  brief  description  of  the  normal  operating  procedure  for  this 
system  follows. 
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Initially,  hydrogen  peroxide  (2  L  of  30%  by  weight  aqueous  solu¬ 
tion)  is  poured  into  the  generator  flask  D  (12  L  spherical)  with  the 
system  open  to  the  atmosphere.  This  solution  is  then  evaporatlvely 
cooled  to  about  15° C  by  evacuating  the  system  to  10-15  torr.  Vapour  is 
solidified  in  the  liquid  nitrogen  trap  H  before  the  vacuum  pump  I. 

A  sodium  hydroxide  solution  (500  n£L,  40%  by  weight)  is  slowly  allowed 
to  enter  the  generator  from  the  storage  vessel  E  (1  L  sphere)  while  a 
minimal  flow  (100  SCCM)  of  inert  gas  (argon)  from  B  enters  the  solution 
through  a  bubbler.  This  bubbler  was  a  20  cm  length  of  flexible  plastic 
(Tygon)  tubing,  sealed  at  one  end  and  pierced  by  about  100  holes  of 
approximately  0.5  mm  diameter.  Trap  C  prevents  liquid  from  flowing 
back,  into  the  gas  handling  system.  Argon  flow  during  the  hydroxide 
addition  ensures  good  mixing.  This  mixing  is  exothermic,  and  the  solu¬ 
tion  temperature  should  be  maintained  below  20°C  (to  minimize  H202 
decomposition)  by  keeping  the  pressure  below  15  torr.  Under  these 
conditions  the  mixing  process  takes  about  20  min. 

02(a)  is  generated  by  flowing  chlorine  through  the  alkaline 
peroxide  solution.  The  reaction  (eq.  1)  is  exothermic.  As  the  solu¬ 
tion  temperature  increases  so  does  the  amount  of  ice  collected  in  the 
liquid  nitrogen  trap.  The  water  vapour  pressure  can,  however,  be  main¬ 
tained  in  the  1-2  torr  range  if,  before  starting  the  C%2  flow,  the 

solution  is  cooled  to  about  -15°C  by  fully  opening  the  system  to  the 
vacuum  pump  for  about  30  min.  Coolant  may  also  be  placed  in  trap  F  of 
Fig.  1  if  it  is  important  to  minimize  H20  vapour  at  the  measurement 
cavity. 

This  system,  operating  with  chlorine  flow  rates  up  to  4000  SCCM, 
produces  total  pressures  up  to  about  10  torr,  of  which  (at  -10*0) 
approximately  1  torr  is  water  vapour.  The  generator  can  operate  con¬ 
tinuously  for  about  40  min  at  1000  SCCM  C&2  and  a  total  pressure  of 

about  3  torr.  After  that,  the  solution  hydroxide  ion  concentration  and 

02(a)  production  simultaneously  fall  rapidly,  and  C%2  begins  to  pass 
through  the  generator. 


>  v. 
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2.2  Problem s  and  Hazards 


The  hazards  of  working  with  C*2  are  well  known  (Refs.  14,  IS). 
The  very  low  recommended  working  threshold  concentration  of  1  ppm  is 
below  the  minimal  (3.5  ppm)  normally  detectable  by  its  odour.  Working 
areas  must  therefore  be  well  ventilated,  and  precautions  must  be  taken 
when  disposing  of  the  generally  small  quantities  of  chlorine  that  con¬ 
dense  in  the  liquid  nitrogen  trap. 

Highly  concentrated  hydrogen  peroxide  also  presents  serious 
safety  hazards,  which  are  well  documented  (Refs.  16,  17).  At  greater 
than  52%  (by  weight)  H202,  the  rate  at  which  heat  is  generated  by 
decomposition  may  exceed  that  at  which  it  is  dissipated  by  evaporation 
of  water  from  the  solution.  The  explosive  generation  of  large  volumes 
of  hot  02/H20  vapour  can  result.  For  this  reason,  30%  H202  has  been 
used  here  rather  than  the  90%  solution  of  the  original  chemical  gener¬ 
ators  (Refs.  2-6). 

Chloride  ion  is  known  to  catalyse  the  decomposition  of  H202. 

The  reaction  sequence  (given  in  Ref.  17,  p.  476)  is: 

H202  +  2C4"  +  2H+  -  C*2  +  2H20  [2a] 

H202  +  Ci2  -  2Ci-  +  2H+  +  02  [2b] 

Therefore,  at  the  end  of  an  02(a)  generation  experiment,  the 
reactor  solution  is  in  an  unstable  condition  since  0H“  produces  Ct~  in 
reaction  1.  As  the  solution  warn,  the  rate  of  the  reaction  sequence 
[2a]-[2b]  increases,  heat  is  evolved,  and  a  violent  decomposition 
ensues.  This  has  occurred  in  our  laboratory  when  the  used  reactor 
solution  has  been  left  overnight  in  the  generator.  It  may  be  avoided 
by  flushing  away  the  still  cool  (<  20* C)  solution  with  a  large  volume 
of  cold  water. 


'."k  .  -  . 
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Removal  of  water  vapour  from  the  gas  flew  is  a  major  problem. 
Even  with  solution  temperatures  in  the  -10®C  to  -15°C  range,  cold  traps 
may  become  blocked  with  ice  after  10-30  min  of  operation,  depending  on 
trap  and  solution  temperatures  and  gas  flow  rates.  At  trap  tempera¬ 
tures  below  about  -70°C  and  gas  temperatures  below  -35°C,  there  is  also 
a  significant  deactivation  of  the  excited  oxygen  (Ref.  6).  A  balance 
must  be  sought  between  permissible  H20  content  and  02(a)  deactivation. 

Degradation  of  the  vacuum  pump  oil  must  also  be  dealt  with.  The 
normal  hydrocarbon-based  pump  oils  are  rapidly  attacked  by  the  acti¬ 
vated  oxygen  (which  is  not  condensed  in  liquid  nitrogen  at  pressures  in 
the  10  torr  range),  and  also  by  traces  of  chlorine  (most  of  which  is 
solidified  in  the  trap).  These  oils  become  viscous  and  datk  brown 
after  about  one  hour  of  continuous  02(a)  generation.  The  risk  of  pump 
oil  explosion  cannot  be  ignored,  particularly  when  other  reactive  gases 
(e.g.  halogens,  nitrogen  fluorides)  may  be  flowing  during  a  laser 
experiment.  Hydrocarbon  oils  should  be  replaced  after  each  such  exper¬ 
iment.  Preferably,  a  completely  inert  pump  fluid,  such  as  Inland  41, 
should  be  used.  For  a  pump  of  several  litres  capacity,  this  can  be 
justified  from  considerations  of  both  safety  and  economy,  since  this 
oil  may  remain  unchanged  during  several  months  of  continuous  usage. 

3.0  MEASUREMENT 

Previously  (Refs.  2-8,  18-22),  electron  paramagnetic  resonance 
(EPR),  isothermal  calorimetry,  and  photometry  have  been  used  to  esti¬ 
mate  absolute  concentrations  of  02(a).  Of  these,  EPR  and  photometry  do 
not  perturb  the  gas  flow  and  are  specific  to  the  singlet  delta  state, 
whereas  isothermal  calorimetry  requires  the  insertion  of  a  probe  into 
the  flow  and  may  be  less  specific  if  other  excited  species  are  pres¬ 
ent. 


,r~.  lie1. 
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We  chose  to  measure  02(a)  using  photon  detection  techniques  for 
the  above  reasons,  as  well  as  because  they  can  follow  concentration 
changes  in  both  space  and  time  during  a  reaction  with  other  species, 
and  are  inexpensive  compared  with  the  EPR  method.  The  disadvantage  is 
that  accuracy  may  not  be  very  high  since  it  depends  on  uncertainties  in 
detector  calibration  and  a  radiative  rate  constant. 

02(a)  emits  most  strongly  in  two  spectral  regions:  in  the  near 
Infrared  around  1268  nm,  and  in  the  visible  centred  at  634  nm  (red). 

The  former  emission  arises  from  direct  spontaneous  transitions  between 
the  vibrational  ground  states  of  the  a  and  X  electronic  states  (the  0-0 
band  of  the  a*4g  -  X3£g  system).  Under  low  resolution,  this  band  has  a 
roughly  Gaussian  distribution  with  a  full  width  at  half  maximum  (FWHM) 
of  about  15  nm  (Ref.  23).  At  higher  resolution  the  rotational  fine 
structure  has  been  partially  resolved  (Ref.  24).  The  red  emission 
results  from  a  collisional  process  involving  two  02(a)  molecules  (hence 
the  common  name  'dlmole'  emission)  that  simultaneously  transfer  their 
electronic  energy  to  one  emitted  photon.  This  dimole  emission  wave¬ 
length  is  thus  half  that  of  the  direct  a-X  emission.  The  dlmole  band 
is  diffuse  with  no  resolved  fine  structure  (Refs.  24-26)  and  an  FWHM  of 
15  nm. 


These  two  emission  processes  will  be  represented  by: 


k 

02(a,0)  +8 

02(X,0)  +  hv(1268) 

[3] 

kd 

202(a,0) 

202(X,0)  +  hv(634) 

[4] 

The  second  bracketed  index  indicates  the  vibrational  level 
(v  ■  0  in  the  cases  indicated)  within  the  electronic  states  a  and  X. 
The  a-state  v'  ■  1  level  is  insignificantly  populated  at  room  tempera¬ 
ture  (assuming  a  Boltzmann  distribution),  but  transitions  to  the  X- 
state  v"  ■  1  level  do  occur: 
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02(a,0)  ♦  02(X,1)  +  hv(1580)  [5] 

202(a,0)  ♦  02(X,1)  +  02(X,0)  +  hv(703)  [6] 

The  second  dimole  emission  at  703  nm,  Indicated  by  [6],  is  of 
comparable  intensity  (Refs.  25,  27)  to  that  at  634  nm,  whereas  the 
1580  nm  emission  of  [5]  is  46  times  lower  in  intensity  than  the  0-0 
band  (Ref.  23). 

Radiative  rate  constants  for  processes  [3]  and  [4]  have  been 
reported  in  the  literature  (Refs.  11,  28  and  18-20,  29  respectively). 
These  constants  directly  relate  the  rate  of  emission  of  photons  in  the 
1268  nm  and  634  nm  bands  to  [o2(a)j  (the  square  brackets  indicate  con¬ 
centration)  : 


dhv(1268)/dt  -  k  [o2(a) ] 

8 

[7] 

dhv(634)/dt  -  kd[02(a)]2 

[8] 

In  [7]  and  [8]  dhv(X„)/dt  is  the  number  of  photons  emitted  per 
second  in  the  band  centred  at  A0,  from  a  volume  of  one  cubic  centi¬ 
meter.  Independent  photometric  methods  based  on  these  emissions  can 
therefore  be  designed  to  estimate  [o2(a)]  without  the  need  to  consider 
other  competing  processes,  such  as  [5]  and  [6],  or  the  presence  of 
other  excited  species  (02(b)  for  example). 

We  have  observed  both  these  emission  bands  with  narrow-band 
interference  filters  and  calibrated  semiconductor  detectors.  The  fol¬ 
lowing  sections  (3.1  and  3.2)  describe  the  apparatus  used  to  produce 
voltage  measurements  from  these  detectors,  and  subsequent  estimates  of 
singlet  delta  concentrations  in  the  02  generator  gas  flow. 
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3.1  Apparatus 

The  measurement  cavity  (G  of  Fig.  1)  of  the  flaw  system  was 
essentially  a  rectangular  plastic  box  (described  in  Section  2.1)  with 
planar  Supracil  II  glass  windows  on  the  sides  parallel  to  the  gas  flow. 
The  viewing  axis  for  the  detectors  and  optics  was  thus  perpendicular  to 
the  flow,  and  since  the  windows  were  sufficiently  long,  the  field  of 
view  was  unconstrained  in  the  flow  direction.  A  separate  report 
(Ref.  30)  presents  a  detailed  development  of  the  equations  that 
describe  the  light  collection  from  this  particular  extended  source,  for 
the  two  quite  different  optical  arrangements  that  we  have  used  experi¬ 
mentally.  Namely,  with  (a)  a  circular  focusing  lens,  and  (b)  a  cylin¬ 
drical  tube  placed  between  source  and  detector.  The  two  arrangements 
provided  independent  methods  for  accurately  defining  the  contributing 
volume  of  emitting  molecules,  and  hence  estimating  their  concentra¬ 
tion. 


Detectors  were  chosen  for  their  high  sensitivity  In  the  two 
emission  regions  of  Interest:  for  the  634  nm  band,  an  EG  &  G  silicon 
HUV-4000B  was  used;  for  1268  nm,  it  was  a  Judson  Infrared  germanium  J- 
16.  Appropriate  narrow-band  interference  filters  (the  transmission 
characteristics  of  which  are  described  in  Appendix  A)  were  placed 
directly  in  front  of  the  detector  surfaces,  with  care  being  taken  to 
exclude  all  stray  light.  A  chopper/lock-ln  amplifier  combination  was 
used  to  record  accurately  the  voltage  signals  from  these  detectors. 

The  detector/filter  units  were  calibrated  absolutely  using  both  an  NBS 
standardized  quartz-iodine  lamp,  and  two  blackbody  sources  at  1000° C. 

3.2  Results 


Table  I  presents  data  from  three  experiments.  For  each  experi¬ 
ment,  column  one  shews  a  series  of  increasing  total  oxygen  pressures 
(P02(total))  measured  at  the  cavity  G  of  Fig.  1.  This  total  02  pres¬ 
sure  was  controlled  by  increasing  the  mass  flow  rate  of  Cl2  through  the 
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generator,  with  the  pump  operating  at  Its  maximum  capacity  (all  pump 
control  valves  fully  open).  The  02  pressure  was  estimated  by  measuring 
the  pressure  in  the  system  Immediately  before  (giving  H20  +  H202  vapour 
pressure)  and  after  the  C&2  flow  was  started,  and  then  again  before  and 
after  the  flow  was  stopped.  The  average  of  these  pressure  differences 
was  taken  to  be  a  measure  of  P02  (total).  This  assumes  that  no  C&2 
passed  through  the  generator  solution,  which  is  supported  by  the  fact 
that  no  significant  C&2  was  collected  in  the  liquid  nitrogen  trap  (H  of 
Fig.  1)  at  the  flow  rates  used  in  these  experiments.  The  estimated  P02 
(total)  was  quite  linear  with  C&2  mass  flow:  e. g.  750  SCCM  C&2  gave 
~  1  torr  02,  1500  SCCM  gave  ~  2  torr,  etc. 

Column  two  of  Table  I  shows  voltages  recorded  at  the  lock-in 
amplifier  from  the  germanium  and  silicon  detector/filter  combinations 
for  the  different  total  oxygen  pressures.  For  the  lens  and  tube  opti¬ 
cal  arrangements  respectively,  eqs.  81  and  82  of  Ref.  30  relate  the 
detector  voltages  from  the  spontaneous  emission  (1268  nm)  to  the  oxygen 
singlet  delta  concentrations,  thus: 

(1)  lens  case  (1268  nm  -  germanium), 

[o2(a) ]  -  4*rv/(wtf)(*tf)(ftf)k  IT  [9] 

8  L 

(ii)  tube  case  (1268  nm  -  germanium), 

|02(a)]  -  4rv/(wtf)(ftf)ksITr2  [lO] 

For  the  dimole  emission  (634  nm),  the  photon  emission  rate 
depends  on  the  square  of  the  02(a)  concentration  (eq.  8),  and  thus  [lo] 
becomes: 

(ill)  tube  case  (634  nm  -  silicon). 
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TABLE  I 


en  singlet  delta  concentrations  from  detector  voltages 


P02 

(torr) 

Voltage 

(Microvolts) 

[02(a)] 

(Molecules,  cm-3) 

P02(a) 

(torr) 

%  0£(b) 

1)  Geraanlus  detector  4  lens: 

1.01 

186 

7.57  El 5 

0.230 

22.8 

2.08 

340 

1.38  E16 

0.421 

20.2 

3.24 

486 

1.98  El 6 

0.601 

18.6 

4.00 

516 

2.10  E16 

0.639 

16.0 

1  1 

11)  Germanium  detector  +  tube: 

1.16 

11.5 

8.28  E15 

0.252 

21.7 

2.19 

18.5 

1.33  E16 

0.405 

18.5 

3.26 

25.5 

1.84  E16 

0.558 

17.1 

3.99 

27.5 

1.98  El  6 

0.602 

15.1 

1  1 

111)  Silicon  detector  +  tube: 

1.02 

0.30 

9.26  El 5 

0.281 

27.6 

2.03 

0.88 

1.59  E16 

0.482 

23.7 

3.29 

1.53 

2.09  El 6 

0.636 

19.3 

4.44 

1.95 

2.36  E16 

0.718 

16.2 

6.00 

2.40 

2.62  El 6 

0.796 

13.3 

A¥ 
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l02(a)]  -  [4rv/(wtf)(ftf)kdITr2]*  [ll] 

The  terms  in  eqs.  9  to  11  are: 

T  -  detector  calibration  constant  (chopping  frequency  23  Hz), 
germanium  (1268  na  filter)  -  1.3E14  photons  s-1v-1; 
silicon  (634  m  filter)  -  9.8E12  photons  s_lv_1; 

V  -  detector  voltage  measured  at  lock-in  amplifier  (in  volts); 

wtf,  itf,  ftf  -  window,  lens,  and  filter  transmission  factors  at  the 
appropriate  wavelength  (see  Appendix  A); 

kg  and  k^  -  spontaneous  and  dimole  emission  rate  constants, 
kg  -  2.6E-4  a-1  (Refs.  11,  28), 

k^  *  4.1E-23  cm3  molecule- *8”*  (an  average  of  three 
reported  values  (Ref.  18)); 

1^  and  IT  -  integrals  over  the  extended  source  volume  for  the  lens  and 
tube  cases  respectively  (their  calculation  is  described  in 
detail  in  Ref.  30); 

r  -  active  detector  radius. 

Equations  9  to  11  lead  to  the  values  for  [o2(a)]  (molecules  cm-3) 
given  in  column  3  of  Table  I.  The  table  also  shows  these  values  in 
torr  (at  20SC,  1  torr  ■  3.29E16  molecules  cm-3)  so  that  the  percentage 
of  02(a)  in  the  measured  total  oxygen  can  be  readily  estimated. 

It  can  be  seen  that  by  increasing  the  total  oxygen  pressure 
(increasing  the  Ci2  flow)  the  absolute  concentration  of  02(a)  arriving 
at  the  cavity  can  be  Increased,  but  that  in  doing  this  the  percentage 
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of  02(a)  in  the  stream  Is  reduced.  This  is  largely  due  to  the  energy¬ 
pooling  process  mentioned  in  Chapter  1.0: 

k 

202(a)  -*-P  02(X)  +  02(b)  [12] 

The  rate  of  this  reaction  will  be  discussed  in  Chapter  4.0 

The  absolute  accuracy  of  the  02(a)  concentrations  can  be  esti¬ 
mated  from  eqs.  9  to  11.  The  detector  calibrations  (T)  and  extended 
source  Integrals  (1^,  1^)  each  have  uncertainties  of  about  10Z,  while 
those  for  the  voltages  and  transmission  factors  are  on  the  order  of  1- 

2Z.  The  error  in  the  radiative  rate  constants  (k  ,  k.)  is  more  dlffl- 

s  d 

cult  to  estimate.  The  value  of  kg  given  by  Badger  et  al  (Ref.  8), 
2.6E-4  s-1,  is  generally  accepted  in  the  literature,  and  agrees  within 
10Z  with  other  quoted  values  (Refs.  8,  28).  The  dimole  emission  rate 
constant  k^  has  been  directly  measured  by  four  independent  groups 
(Ref 8.  18-20,  29).  Of  these,  three  agree  to  within  10Z  (Ref.  18), 


whereas  the  value  found  by  Derwent  and  Thrush  (Ref.  20)  is  smaller  by  a 
factor  of  two. 

An  uncertainty  of  about  30Z  is  therefore  expected  in  the  [o2(a)] 
values  given  in  Table  I  from  observation  of  the  direct  spontaneous 
emission  (germanium  detector). 


vmjn 


The  concentrations  estimated  from  the  dlmole  emission  (silicon 
detector)  have  an  Inherently  higher  accuracy  due  to  the  square  root 
relationship  in  [ll]«  The  values  from  the  silicon  detector  in  Table  I 
therefore  have  an  uncertainty  of  about  15Z. 


Figure  2  shows  the  calculated  partial  pressures  of  02(a) 
arriving  at  the  measurement  cavity  for  each  of  the  measured  total  oxy¬ 
gen  pressures  in  the  two  detector  data  sets. 


/  ,*  ** 
i  **• 


-  *  ■.a*j**« 


vV.v 

VmVmV, 


^ .  •  -  - .  •  •  ■  -  • ,  >  •  -  •  -  *  e  v-  v  v  v  v  v*  VvV;,*  >  v  *.*  *.* 
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□ 


Q  Silicon  Detector  Data 
O  Germanium  Detector  Data 
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It  Is  significant  that  the  [o2(a)]  values  calculated  from  the 
germanium  detector  data  were  consistently  lower  than  those  calculated 
from  the  silicon  detector  voltages.  The  differences  are  somewhat 
greater  than  would  be  expected  from  uncertainties  in  the  detector  cali¬ 
brations  alone,  and  suggest  that  the  accepted  value  for  the  spontaneous 
radiative  rate  constant  (k^  ■  2.6E-4  s”1)  for  02(a)  may  be  rather  high. 
The  dimole  emission  rate  constant  has  been  more  extensively  studied  and 
is  probably  more  reliable. 

The  curve  drawn  through  the  experimental  points  has  the  func¬ 
tional  form: 


y(x)  ■  x/(a  +  bx)  [ 13 ] 

(y  =  P02(a);  x  =  P02(total)). 

This  form  was  chosen  because  it  has  the  correct  asymptotic 
behaviour  (y  *  constant  (1/b)  as  x  ♦  *)  for  a  second-order-dominated 
decay  process,  and  it  leads  to  a  better  least-squares  fit  than  the 
other  simple  two-parameter  forms  axb  and  x(a  +  bx).  Since  there  are 
eight  data  points  from  the  germanium  detector  and  only  five  from  the 
silicon,  a  weighted  fit  was  performed  giving  the  silicon  points  1.6 
times  the  weight  of  the  germanium  points.  This  means  that  the  data 
sets  from  each  detector  were  treated  equally. 

The  resulting  values  for  the  constants  of  [ 13 ]  are: 

a  -  3.213  ;  b  -  0.7449  [l4] 

No  physical  significance  is  to  be  attached  to  these  values. 
Equations  13  and  14  merely  provide  a  continuous  representation  of  the 
02(a)  partial  pressure  at  the  cavity  as  a  function  of  total  02  pres- 


UNCLASSIFIED 

17 


3.3  Inert  Gas  Diluents 


We  have  studied  the  effect  on  the  02(a)  emissions  of  introducing 
up  to  10  torr  of  either  argon  or  sulphur  hexafluoride  into  the  gas 
stream.  The  diluent  gas  was  premixed  with  the  C and  passed  directly 
through  the  generator  solution. 

No  detectable  effect  has  been  observed  on  either  the  1268  ran  or 
634  raa  emission  Intensities  for  any  fixed  CJt2  flow  rate.  These  obser¬ 
vations  are  in  accord  with  the  reported  (Refs.  31,  32)  low  02(a) 
quenching  rate  constants  (IE-20  cm3  molecule-1  s-1)  for  these  gases. 
Inert  behaviour  can  also  be  expected  for  nitrogen  and  helium  on  the 
basis  of  their  similarly  low  quenching  rates. 

3.4  Water  Vapour 

At  a  generator  solution  temperature  of  -10°C  there  was  about 
1  torr  of  solution  vapour  in  the  flow  system,  measured  at  the  cavity. 
This  was  essentially  all  water  vapour  Blnce  the  mole  fraction  of  H202 
in  the  vapour  above  a  30%  aqueous  solution  is  <  0.003  (Ref.  17,  page 
227)  at  this  temperature.  As  Ci2  flows  into  the  reactor  the  solution 
temperature  rises  slowly  (about  1*C  per  minute  for  2000  SCCM  of  Ci2), 
and  the  measured  emission  intensity  declines  steadily  as  this  tempera¬ 
ture  exceeds  approximately  0°C. 

The  02(a)  quenching  rate  constant  for  is  about  5E-18  cm-3 
molecule-1  s- 1  (Refs.  31,  32),  which  is  sufficiently  high  to  produce  a 
measurable  effect  with  increasing  H20  vapour  pressure  (see 
Chapter  4.0).  Operating  the  generator  at  a  solution  temperature 
between  -10°C  and  -5°C  (in  the  1-2  torr  H20  pressure  range)  does  not 
give  a  significant  variation  in  emission  Intensity  for  a  fixed  C&2 
flow.  The  measurements  given  in  Table  I  were  all  taken  with  the  solu¬ 
tion  in  this  temperature  range. 
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3.5  Vapour  Removal  by  Liquid  Nitrogen  Trap 

The  effect  of  trap  temperature  on  the  percentage  02(a)  In  the 
flow  has  been  studied  extensively  elsewhere  (Refs.  6,  7).  We  used 
liquid  nitrogen  in  trap  F  of  Fig.  1  to  study  how  almost  entirely  remov¬ 
ing  H20  from  the  flow  affects  reactions  with  NF  radicals.  Under  these 
conditions  the  percentage  of  02(a)  is  dramatically  reduced  to  between  3 
and  5Z  for  the  pressure  range  1-4  torr  total  02  (cf.  25-16Z  in  Fig.  2 
for  these  pressures).  Trap  temperatures  as  low  as  -78°C  (acetone/dry 
ice)  are  possible  (Ref.  7)  without  significant  quenching  of  02(a). 

4.0  DEACTIVATION 


First  and  second  order  rate  constants  quoted  in  this  and  subse¬ 
quent  sections  have  the  units  s“*  and  cm3  molecule-*  s-1  respectively. 
Concentrations  are  in  molecules  cm-3  unless  otherwise  stated. 

There  are  three  principal  mechanisms  for  the  deactivation  of 

02(a) . 


(i)  Energy  Pooling: 

2A  B  +  X 

[15] 

(ii)  Quenching: 

kA 

A  +  Qi  ™  X  +  Qi 

[16] 

(ill)  Wall  Collisions: 

kA 

A  X 

[17] 

In  [ 15 ]— [ 17 ]  A,  B  and  X  represent  02(a),  02(b)  and  02(X)  respec¬ 
tively,  and  Q1  is  any  other  gaseous  molecule  that  may  be  present  (H20, 
Ar,  and  02(x)).  The  same  alphabetic  symbol  will  represent  species 
concentration  whenever  there  is  no  ambiguity. 
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Radiative  losses  (eqs.  3  to  6)  are  negligible  in  the  time  domain 
considered  here  (<  Is). 


In  this  section,  an  equation  determining  the  rate  of  change  of 
[o2(a>]  Is  developed,  and  the  Integrated  form  is  derived.  The  02(a) 
concentration  Is  then  defined  for  all  time  in  terms  of  the  initial  gas 
concentrations  and  the  reaction  rate  constants.  The  [o2(a)j  curve  of 
Fig.  2  can  therefore  be  reproduced,  given  the  flow  time  from  generator 
to  measurement  cavity,  using  a  particular  set  of  rate  constants. 

The  rate  of  loss  of  A  due  to  reactions  [  15 ]— [  17 ]  is: 

-A'(t)  -  2kpA2  +  (k*  +  J  k*£  Qi)A  [18] 

Note  that  the  sun  over  1  includes  all  quenching  species  in  the  flow, 

and  that  the  rate  of  loss  of  A  due  to  pooling  is  written  2k  A2  in 

P 

agreement  with  the  literature  (Refs.  20,  21,  36,  37)  definition  of  kp 
(which  is  the  rate  constant  for  the  formation  of  B  and  X). 


There  is  a  great  deal  of  evidence  in  the  literature  (Refs.  37- 
44)  to  indicate  that  the  02(b)  produced  by  the  pooling  reaction  [ 15 J  is 
quenched  to  02(a)  both  by  other  gaseous  species,  and  wall  collisions. 
Reactions  that  produce  A  must  therefore  be  included: 


A 


[19] 


[21] 
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where: 

14  -  *£  +  £  kqjQi  5  (I  -  A,  B)  [22] 

The  rate  equation  for  B  due  to  reactions  [ 15 ] »  [ 19 ]  and  [20] 
la: 

B‘(t)  -  kpA2  -  KbB  [23] 

Water  vapour  Is  abundant  In  the  chemical  generation  of  02(a), 

and  It  la  an  extremely  rapid  quencher  of  02(b).  Indeed,  Q  is  about 

5E-12  (Refs.  31,  38,  45,  46),  so  that  for  1.5  torr  H20  the  pseudo-first 

order' quenching  rate  la  roughly  250,000  s-1 .  Other  quenchers,  and  wall 

collisions,  are  completely  negligible  In  the  term  for  14  (eqs.  22  and 

23):  is,  for  example,  approximately  4E-17  (Refs.  38,  47-49). 

u2 

02(a),  however,  la  quenched  about  a  million  times  more  slowly  by  water 
vapour  than  la  02(b)  (Refs.  31,  32).  This  rapid  B  quenching  appears  to 
be  due  to  a  near  resonance  between  the  02  b-a  energy  spacing  and  a 
vibrational  mode  of  H20. 


The  effect  of  this  very  fast  conversion  of  B  to  A  is  that  a 
quas 1-steady-state  for  B  Is  established  very  rapidly:  for  1  torr  of 
Initial  A,  B  rises  to  an  essentially  constant  value  in  ~  5E-5  s.  This 
was  calculated  by  exactly  solving  the  coupled  eqs.  21  and  23  with  rate 
constants  from  the  literature  using  the  general  chemical  kinetics 
FORTRAN  program  described  In  Ref.  50.  The  quasl-steady-state  concen¬ 
tration  of  B  is  obtained  from  [23]: 
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A* (t)  -  “kpA2  “  KaA  [25] 

A(t)  is  still  coupled  to  X(t)  through  the  term  in  KA  due  to 
the  quenching  of  A  by  X  (cf.  eq.  16).  The  case  Qi  -  X  will  be  sepa¬ 
rated  from  the  general  quenching  summation: 

-A'(t)  -  kpA2  +  (k*  +  k*X  ♦  J  k*±  Qi)A  [26] 
and  X(t)  will  be  eliminated  by  assuming 

X(t)  -  X(0)  +  A(0)  -  A(t)  =  (X  +  A  )  -  A(t)  [27] 

00 

Equation  27  is  valid  because  B(t)  is  always  negligible  in  a  flow 
containing  more  than  1  torr  of  water  vapour  (cf.  [24]).  XQ,  the 
Initial  ground  state  O2  concentration,  is  essentially  zero  for  the 
chemical  generation  process  since  there  is  close  to  100%  02(a)  at  the 
llquid/gas  Interface  in  the  generator  (Ref.  8).  Deviations  from  this 
condition  will  be  investigated. 

When  [27]  is  used  for  X  in  [26]  a  completely  uncoupled  first- 
order  differential  equation  for  A(t)  results.  This  may  be  Integrated 
directly  to  give: 


where: 


A(t> 


A  K_ 
o  T 


(KT  +  Vo)  exp  <V>  "  Vc 


[28] 


kD  =  kp  “  '‘x 


s 5  k» +  >4  <*« +  V  ^  kQi  01 


[29] 

[30] 


LA* 
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The  superscript  A  on  the  rate  constants  will  be  omitted  In  sub¬ 
sequent  discussion  since  all  constants  involving  B  have  been  elimi¬ 
nated* 


Estimation  of  the  Pooling  Rate  Constant 


Although  A  has  been  written  A(t)  in  [28],  it  is  more  generally 
the  many-variable  function  A(t,  Aq,  Xq,  {k^,  Qi}»  kw>  k  ).  So  for  a 
fixed  time  (at  a  fixed  region  in  space  in  the  flow  system:  the  meas¬ 
urement  cavity),  [28]  can  generate  02(a)  concentration  as  a  function  of 
total  o^rgen,  given  a  set  of  values  for  the  rate  constants  and 
quenching  species  concentrations.  In  particular,  the  experimental 
curve  of  Fig.  2  can  be  matched  with  a  theoretical  curve  using  values 


for  kw,  kp,  kx,  k^  and 
system) . 


(partial  pressure  of  water  vapour  in  the 


An  estimate  of  the  flow  time  between  generator  and  cavity  is 
required.  We  made  a  direct  measurement  of  this  time  using  two 
detector/chopper/lock-ln  amplifier  combinations.  One  detector  viewed 
the  liquid-gas  interface  in  the  generator,  while  the  other  was  at  the 
downstream  measurement  cavity.  Their  voltage  outputs  were  recorded  on 
a  dual-pen  chart  recorder,  whose  paper  was  set  to  move  at  40  in/min. 

On  opening  the  chlorine  flow,  we  recorded  a  time  delay  between  the 
detector  signals  of  0.43  s,  reproducible  to  within  0.06  s.  We  there¬ 
fore  take: 
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The  wall  deactivation  constant  kw  was  found  to  vary  considerably 
(Refs.  31-35)  depending  on  the  condition  and  material  of  the  walls.  A 
glass  tube  freshly  rinsed  with  HF,  and  evacuated  for  several  hours,  may 
have  kw  as  low  as  0.1  (Ref.  35),  whereas  others  have  reported  values  In 
the  range  0.2-0. 6  (Refs.  32,  8).  Only  at  higher  pressures  (above 
100  torr)  does  kw  become  pressure  dependent  (Ref.  32)  (diffusion 
limited).  The  wall  deactivation  constant  has  therefore  been  treated  In 
this  study  as  a  variable  parameter.  Since  during  the  chemical 
generation  process  the  walls  near  the  generator  become  coated  with  a 
fine  solid  deposit  (NaOH,  NaCi),  a  higher  wall  deactivation  rate  can  be 
expected. 

The  generally  accepted  literature  value  for  the  pooling  rate 
constant  is  that  of  Derwent  and  Thrush  (Ref.  20),  who  gave 
kp  ■  (2.0  ±  0.5)E-17.  Recently,  Fisk  and  Hays  (Ref.  37)  also  reported 
a  value  of  2.0E-17  for  k^,  and  gave  an  estimated  uncertainty  of  ±  302. 

A  value  of  2.4E-17  has  also  been  quoted  (Ref.  53).  However,  two 
Independent  groups  (Refs.  21,  36)  had  previously  obtained  values  about 
ten  times  lower  than  this  (0.22,  0.23E-17,  respectively).  The  first  of 
these  groups  (Arnold  and  Ogryzlo)  had  also  estimated  an  upper  limit  of 
1.0E-17 . 

In  this  study,  fixed  values  of  k  between  0.2E-17  and  3.0E-17 
were  investigated.  For  each  value  of  k^,  eq.  28  was  used  to  generate  a 
numerical  least-squares  fit  to  points  on  the  experimental  curve  of 
Fig.  2  by  variation  of  the  constant  k^.  Twelve  points  were  chosen  at 
0.5  torr  Intervals  In  the  total  02  pressure  range  0.5  to  6  torr.  Thus, 
for  a  given  k^  an  optimum  curve  can  be  obtained  by  numerical  minimiza¬ 
tion  (wrt  kw)  of  the  least-squares  function: 


12 r  expt.  theory  .. 

S  =iIi[(P02(a))1  -  (P02(a))1]2 


[37] 
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A  minor  problem  arises  when  eq.  28  is  used  with  the  experimental  total 
02  pressures  of  Fig.  2.  These  were  pressures  ^  the  measurement  cav¬ 
ity,  whereas  the  pressures  In  the  generator/flow  system  up  to  the 
cavity  were  consistently  6X  higher  because  there  was  a  short  (10  cm) 
length  of  smaller  diameter  tubing  (2.5  cm)  connecting  the  two  regions. 
The  effect  of  this  small  pressure  difference  has  been  Included  in  the 
FORTRAN  program  that  performs  the  numerical  least-squares  fits. 

The  form  of  the  minimized  S  Is  shown  In  Fig.  3  over  the  range  of 

values  of  k  .  There  Is  a  well-defined  minimum  at 
P 

k  -  2.23E-17  [38] 

P  1 

for  which  the  optimum  wall  deactivation  rate  constant  Is 

k  -  2.55  [39] 

w  1  J 

This  pair  of  values  (k  ,  k  )  therefore  gives  the  best  possible 

p  w 

fit  to  the  experimental  curve  of  Fig.  2. 

Figure  4  shows  our  best-fit  curve  through  the  experimental 
points  taken  from  the  curve  of  Fig.  2.  It  also  shows  the  best  fit 
obtained  using  the  k^  given  by  Arnold  &  Ogryzlo.  It  Is  evident  that 
their  value  of  0.22E-17  for  k^  cannot  reproduce  the  curvature  that  we 
obtained  experimentally. 

4.2  Analysis  of  Errors  In  the  Estimation  of  k 

The  value  of  k^  that  minimizes  S  (eq.  37)  depends  on  the  values 
assigned  to  the  following  parameters: 

1)  the  set  of  oxygen  singlet  delta  concentrations  at  the 

measurement  cavity;  l.e.  the  experimental  values  (po2(a)}  of 
Figs.  2  and  4; 
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2)  the  amount  of  oxygen  that  is  In  the  ground  state  at  the  gen¬ 
erator  solution/gas  interface;  i.e.  Xq  of  eqs.  27  to  30; 

3)  the  oxygen  quenching  rate  constant  k^;  and 

4)  the  time  t  that  molecules  spend  in  the  flew  system  from  gen¬ 
erator  solution  to  measurement  cavity. 

The  other  parameters  required  in  eqs.  28  to  30,  which  are 

quenching  species  concentrations  and  rate  constants,  only  contribute  to 

the  pseudo-first-order  rate  constant  K^.  They  do  not  Influence  the 

computed  value  of  the  second-order  constant  k^.  However,  they  do 

affect  the  optimized  value  of  k^,  which  is  therefore  less  accurately 

known  than  k^.  In  this  study,  only  vapour  has  been  Included  as  a 

quencher.  The  optimized  values  for  k^  therefore  contain  contributions 

from  pseudo-first-order  quenching  by  H202  and  Ci2.  Since  these  have 

low  concentrations  in  the  gas  stream,  the  contributions  are  expected  to 

be  small.  First-order  quenching  due  to  solution  droplets  entrained  in 

the  gas  stream  also  contributes  to  our  value  of  k  . 

w 

Table  II  shows  the  effect  of  the  uncertainties  in  the  input 

parameters  on  the  optimized  values  for  k  and  k  (Ak  ,  Ak  )  when  a 

p  w  p  w 

given  parameter  (column  1)  takes  either  the  maximun  (+)  or  minimum  (-) 
value  defined  by  the  error  limits  given  in  column  3.  Thus,  row  1  shows 
that  when  the  experimentally  measured  02(a)  concentrations  are 
decreased  by  30%  the  best-fit  k^  value  increases  by  0.0296E-17,  a 
change  of  only  13Z.  The  optimized  k^,  on  the  other  hand,  Increases  by 
33Z.  kw  is  much  more  sensitive  to  changes  in  the  absolute  concentra¬ 
tions  of  02(a)  than  is  k^.  It  is  the  curvature  of  P02(a)  as  a  function 
of  increasing  total  oxygen  pressure  that  is  a  sensitive  measure  of  k^ 
(cf.  Fig.  4).  Row  2  indicates  the  effect  of  a  deviation  from  100Z 
production  of  02(a)  at  the  generator.  Initially,  Xq  was  taken  to  be 
zero.  When  10Z  of  the  total  02  at  the  generator  is  taken  to  be  in  the 
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TABLE  II 

Effect  of  Input  parameter  uncertainties  on  kp  and  ky 


Parameter 

Value 

Uncertainty 

4k  x 
P 

10*7 

4k 

w 

(+) 

(-) 

(+) 

(") 

{P02(s) } 

Fig.  4 

±  30Z 

-0.206 

0.296 

-0.612 

0.833 

X. 

0 

+  10Z  02 

0.137 

— 

-0.246 

— 

k 

1.5E-18 

±0.1E-18 

-0.025 

0.025 

0.005 

-0.005 

X 

t 

0.43 

±  0.06 

-0.319 

0.422 

-0.333 

0.441 

^0 

4.8E-18 

±0.8E-18 

— 

— 

-0.039 

0.039 

[h2o] 

1.5  torr 

±0.3  torr 

-0.048 

0.048 

ground  state,  k  Increases  by  6Z,  and  k  by  10Z.  It  has  been  estimated 
P  ^ 

(Ref.  8)  that  at  least  95Z  of  the  oxygen  is  in  the  a-state  at  the 
generator,  so  that  row  2  gives  a  generous  estimate  of  this  source  of 
error. 


Error  limits  for  k  and  k  can  be  established  by  combining  the 
P  ^ 

effects  of  the  uncertainties  in  the  parameters  of  Table  IT.  We  obtain: 


kp  -  2.2  (+  0.9  or  -0.5)  E-17  [40] 


k  -  2.6  ±  1.3 
v 


[41] 


maaamaaBaaams^m 
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The  temperature  of  the  gas  also  has  an  effect.  The  generator 

solution  Is  in  the  -10  to  -5°C  range,  and  the  gas  arrives  at  the  cavity 

at  room  temperature  (20°C).  However,  recent  publications  (Refs.  51, 

52)  show  the  variation  of  k  to  be  small  from  -15  to  20°C.  If  we 

P 

assume  the  average  gas  temperature  to  be  about  0°C,  kp  would  be  reduced 
by  approximately  10%  from  Its  value  at  295  K. 

Our  results  are  therefore  consistent  with  the  recent  values  for 
kp  of  about  2E-17  (Refs.  20,  37,  53). 

5.0  CONCLUSIONS 


Oxygen  singlet  delta  can  be  chemically  generated  in  a  vacuum 
flow  system  without  serious  difficulty.  Safe  operating  techniques  and 
photometric  methods  for  02(a)  concentration  were  developed  for  the 
generator  that  we  constructed  at  DREV . 

Operating  at  a  pressure  of  1  torr  total  oxygen,  this  generator 
produced  (25  ±  7)X  02(a)  at  a  downstream  measurement  cavity.  The  per¬ 
centage  02(a)  decreased  to  about  13%  as  the  total  02  pressure  was 
increased  to  6  torr.  Argon  and  SFfi  used  as  diluent  gases,  up  to  a 
partial  pressure  of  10  torr,  did  not  measurably  alter  this  perform¬ 
ance. 


The  excited  oxygen  Is  deactivated  principally  by  water  vapour, 
wall  collisions,  and  energy-pooling.  It  should  be  possible  to  obtain 
close  to  50%  02(a)  In  a  downstream  cavity  by: 

1)  maintaining  the  generator  solution  temperature  at  -15eC; 

2)  minimising  the  flow-system  volume;  and 

3)  doubling  the  pumping  rate  to  about  20  L/s. 
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We  estimate  the  pooling  rate  constant  for  oxygen  singlet  delta 

to  be: 

k  *  2.2E-17  cm3  molecule-1  s-1  [43] 

P 

with  upper  and  lower  limits  of  3.1E-17  and  1.7E-17  respectively*  For 
the  chemical  generation  process,  a  pseudo-first-order  decay  rate  con¬ 
stant 

k  -  2.6  ±  1.3  s-1  [44] 

w 

is  appropriate. 

With  these  two  values,  eq.  28  may  be  used  to  predict  the  per¬ 
formance  of  any  chemical  oxygen  generator  if  a  reliable  estimate  of  the 
generator  to  measurement  cavity  flew  time  is  available. 

This  work  on  the  chemical  generation  of  oxygen  singlet  delta  has 
given  us  the  potential  to  rapidly  develop  a  chemical  oxygen  iodine 
laser  operating  at  1.315  Mm,  and  to  study  the  blue-green  emission  at 
529  nm  arising  from  a  pooling  reaction  between  the  singlet  delta  states 
of  oxygen  and  nitrogen  fluoride. 
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APPENDIX  A 

Transmission  Factors 

This  section  discusses  the  fraction  of  light,  emitted  In  a  band 
associated  with  some  set  of  molecular  transitions,  that  passes  through 
a  medlm  for  which  the  transmittance  Is  known  as  a  function  of  wave¬ 
length. 

Let  h(A)  be  a  function  of  wavelength  describing  the  relative 
Intensity  of  the  light  emitted  in  such  a  band,  and  let  h(A)6A  be  the 
fraction  of  all  the  photons  In  the  band  that  are  emitted  In  the  Infini¬ 
tesimal  wavelength  interval  $A.  Thus: 

7  h <A)dX  -  1  [A-l] 

o 

If  g(A)  gives  the  transmittance  of  the  medium,  then  the  fraction 
of  photons  from  the  hand  that  passes  through  the  medium  in  the  Interval 
«A  Is  g(A)  h(A)  6A. 

A  transmission  factor,  tf,  may  thus  be  defined  for  the  entire 

band: 

7  g(A)  h(A)  dA  =  tf  [A-2] 

o 

An  emission  profile  under  low  resolution  may  often  be  approxi¬ 
mated  by  a  Gaussian  function,  thus: 

h(A)  -  ho  exp(Y(X-Ao)2)  [A— 3  J 

and  [A-l]  may  be  used  to  define  hQ,  giving 

[A-4] 


h(A)  -  (Y/»)%  exp(-Y(A-A)2) 
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Y  is  given  by  Che  full  width  at  half  maxlmim  (FWHM)  of  the  emis¬ 
sion  peak: 


_  -41n(0. 5) 
(FWHM) z 


[A-5] 


If  the  transmission  medium  is  a  window  or  a  lens,  for  which  the 
transmittance  is  generally  a  constant  (T)  over  the  wavelength  region  of 
the  emission  band,  eq.  A-2  is  trivial: 


(itf)  or  (wtf)  -  / Th(X)dX  -  T  [a-6] 


The  transmittance  of  a  narrow-band  Interference  filter,  on  the 
other  hand,  may  often  be  approximated  by  a  Gaussian: 

g(A)  -  a  exp(-3(X-Xf)2)  [A-7 ] 

where  a  is  the  maximum  transmittance,  at  X  *  Xf,  and  3  is  calculated 
from  the  FWHM  of  the  filter  transmittance  curve  (cf.  eq.  A-5).  Gener¬ 
ally,  the  centre  wavelength  X^  is  chosen  equal  to  X,  for  the  emission 
band.  In  this  case  [a— 2 ] ,  [A-3]  and  [a-7]  give: 

(ftf)  -  a(Y/ir)*  /  e"(Y+0)(X_Xo)2  dX  [a-8] 
o 


l.e. 


(ftf)  -  a(Y/Y+B)*  [A-9] 

In  general,  the  transmittance  curve  for  the  filter  and  the 
emission  band  profile  may  be  treated  polntwise  in  X.  The  integral  in 
[A-2]  can  then  be  calculated  numerically  to  a  high  level  of  accuracy  if 
many  points  are  used. 


UNCLASSIFIED 

37 


In  some  cases  the  filter  transmittance  curve  may  be  well 
approximated  by  a  rectangle  or  a  trapezium. 

In  the  estimates  of  02(a)  concentrations  described  in  this 
report,  two  emission  bands  were  observed  using  different  narrow-band 
Interference  filters. 

(1)  Dimole  emission  at  10  ■  634  nm 

This  was  observed  through  a  filter  centred  at  634.6  nm,  with 
FWHM  *  10.8  nm,  and  a  transmittance  curve  well  approximated  by  a  trape¬ 
zium  (a  **  0.61).  The  emission  band  is  smooth  with  no  resolved  struc¬ 
ture,  FWHM  ■  15  nm  (Refs.  24-26),  and  is  well  represented  by  a  Gaussian 
form  in  the  region  of  significant  filter  transmittance.  The  Integral 
in  [a— 2 ]  was  performed  numerically  to  yield: 

(ftf)634  -  0.364  ±  0.005  [a-10] 

(2)  Emission  at  A,  ■  1268  nm 

A  filter  centred  at  1268  nm,  with  FWHM  -  20  nm  and  a  -  0.54,  was 
used.  A  Gaussian  form  approximated  the  transmittance  curve.  The  emis¬ 
sion  band  at  low  resolution  (Refs.  23,  24)  is  also  approximately 
Gaussian,  with  FWHM  -  15  nm.  Equation  A-9  can  then  be  used,  with  Y  and 
6  given  by  [ A— 5 ] : 

(ftf)126g  -  0.43  ±  0.01  [a-11] 
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